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ABSTRACT 

Wc have analyzed the spatially integrated kinematic properties of the ionized gas 
within the inner r <few kpc in 13 optically selected SDSS type 2 quasars at z ^0.3-0.6, 
using the [OIII] AA4959,5007 lines. The line profiles show a significant asymmetry in 
11 objects. There is a clear preference for blue asymmetries, which are found in 9/13 
quasars at 10% intensity level. In coherence with studies on other types of active and 
non active galaxies, we propose that the asymmetries are produced by outflows where 
differential dust extinction is at work. 

This scenario is favoured by other results we find: in addition to quiescent ambient 
gas, whose kinematic properties are consistent with gravitational motions, we have 
discovered highly perturbed gas in all objects. This gas emits very broad lines (R = 
F p\^H M ° t 111 ^ ^2). While the quiescent gas shows small or null velocity shifts relative 
to the systemic velocity, the highly perturbed gas trends to show larger shifts which, 
moreover, are blueshifts in general. Within a given object, the most perturbed gas 
trends to have the largest blueshift as well. All together support that the perturbed 
gas, which is responsible for the blue asymmetries of the line profiles, is outflowing. 
The outflowing gas is located within the quasar ionization cones, in the narrow line 
region. 

The relative contribution of the outflowing gas to the total [OIII] line flux varies 
from object to object in the range ^10-70%. An anticorrelation is found such that, 
the more perturbed the outflowing gas is, the lower its relative contribution is to the 
total [OIII] flux . This suggests that outflows with more perturbed kinematics involve 
a smaller fraction of the total mass of ionized gas. 

Although some bias affects the sample, we argue that ionized gas outflows are a 
common phenomenon in optically selected type 2 quasars at 0.3< z <0.6. 
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1 INTRODUCTION 

Evidence for an intimate connection between supermassive 
black hole (SMBH) growth and evolution of galaxies is nowa- 
days compelling. Not only have SMBHs been found in many 
galaxies with a bulge component, but tight correlations ex- 
ist between the black hole mass and some bulge properties, 
such as t he ste llar mass and vel ocity dispersion a (Ferrarese 
& Merrit |2000| , Gebhardt et al- bOOd ). The origin of this re- 
lation is still an open question, but quasar induced outflows 
might play a critical role. Hydrodinamycal simulations show 
that the energy input from quasars can regulate the growth 



* Based on observations carried out at the European Southern 
Observatory (Paranal, Chile) with FORS2 on VLT-UT1 (pro- 
grammes 081.B-0129 and 083.B-0381. 



and activity of black holes and their host galaxies (di Mat- 
teo, Springel & Hernquist 120051 ). Such models show that a 
merger between two galaxies with SMBH leads to strong 
inflow that feeds gas to the SMBH activating the quasar 
phase. The energy released by the strong outflows associ- 
ated with major phases of accretion expels enough gas to 
quench both star formation and further black hole growth. 
This determines the lifetime of the quasar phase and ex- 
plains the relationship between the black hole mass and the 
stellar velocity dispersion. 

This scenario, although vey attractive, still lacks com- 
pelling observational support. The presence of outflows in 
different classes of active galaxies has been known for sev- 
eral decades. However, it is still not clear whether they can 
generate important feedback effects with a notable impact 
on the systems evolution. The limitations are several, such 
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as the uncertainty on where the outflows originate and their 
geometry, or the fact that several gas phases (molecular, 
neutral, ionized) can be involved. These are often observable 
at different spectral ranges and spatial scales and therefore 
require different observing techniques. 

The presence of outflows of ionized gas in different types 
of active galaxies is well stablished. They are evident, for 
instance, as absorption lines in the UV and X-ray spec- 
tra, blueshifted relative to the systemic velocity of the host 
galaxy. Outflowing velocities of up to —2000 km s" 1 h ave 
been measured (see Crenshaw, Kraemer & George 1999 for 
a review). 

Outflows of ionized gas have also been discovered when 
characterizing the gas kinematics using the spectral proper- 
ties of a diver sity of forbidden lines in the infared (e.g. Spoon 
& Holt l20ld ) and the optical. The [OIII]AA4959,5007 lines 
have been used more frequently in optical studies. These 
are often the strongest optical lines and have the advantage 
of being cleanly isolated from other emission and absorp- 
tion features. In addition, underlying contamination by line 
emission from the broad line region (BLR) is not expected, 
as is often the case for permitted lines such as Ha and H/3. 

The predominance of blue assymetries on the [OIII] line 
profiles, their frequent blueshift relative to different indica- 
tors of the systemic redshift and the kinematic substructure 
of the lines indicate that the narrow line region (NLR) gas 
is not a stable, quiescent gas reservoir that follows the bulge 
gravitational potential. This dynamic behaviour has often 
been interpreted in terms of outflows where differential red- 
dening is present: i.e. at least part of the NLR gas is moving 
radially outwards and dust reddening is more severe for the 
receding gas. This applies to active galaxies in general (e.g. 
Greene & Ho l2005l) and specific AGN t ypes: Seyferts 1 and 
2s (Whi ttle et al. Il988t Colb ert et al. 1 19961 . Crenshaw et 
al. I2010L Heckman et al. Il98ll . Veille ux ll99ll ). Narro w Lin e 
Seyfert 1 galaxies (Bian et al. 120051 . Komossa et al. 2008), 
radi o loud and radio quie t type 1 quasar s (Le ipski & Ben- 
nert 120061 . Heckman et al. Il984 Boroson 120051 ) and narrow 
line r adio galaxies at differe nt redshifts (Humphrey et al. 
2006 1. Heckman et al. Il98ll. Holt, Tadhunter & Morganti 
2008h . iHeckman et al.l l|l98ll ) found that blue assymmetries 
are in general not found in steep spectrum (extended) ra- 
dio sources, but preferentially in radio quiet objects and flat 
spectrum (compact) radio sources. 

The origin of the outflows has also been subject of study. 
Correlations between the radio properties and the kinematic 
signatures of the outflows have suggested that the radio 
jets/lobes can play an important role to trigger such outflows 
both in rad io loud and radio q uiet o bjects (e.g. Hum phrey 
et al. 120061 . L eipski & Bennert 120061 . Heckman et al. 1 1984 
Whittle et al. Il988l . etc.). The dynamics of the ionized gas, 
which is subject to a partial or full radiation field from the 
central continuum source, may be affected by the radiation 
force as well (e.g. Boroson [2005). Finally, galactic super- 
winds may also exist. Star forming systems at low and high 
redshift are known to drive galactic winds with observable 
sig natur es in the wings of the emission lines (e.g. Shapiro et 
al.lH). 

Type 2 (obscured) quasars have been discovered in large 
quantities only in rece nt yea rs (e.g. Zakamska et al. 120031 . 
Martfnez-Sansigre et al. l2005l ). From the observational point 
of view, these are the luminous counterparts of Seyfert 2 



galaxies. In particular, Zakamska et al. identified in 2003 
several hundred objects in the redshift range 0.3^ z ^CO.8 
in the Sloan Digital Sky Survey (SDSS) with the high ion- 
ization narrow emission line spectra characteristic of type 2 
AGNs and narrow line luminosities typical of type 1 quasars. 
Their far-infrared luminosities place them among the most 
luminous quasars at similar z. Several studies suggest that 
th ey hav e high star-f ormin g luminosities (e.g. Zakamska et 
al. 120081 . Hiner et al. 120091 ). They show a wide range of X- 
ray lu minosities and obscur ing column densities (Ptak et al. 
2006| Zakamska et al. l2004 ). The host galaxies are most fre- 
qu ently ellipticals with irregular morphologies (Zakamska et 
al. 120061 ). and the nuclear optical e missi on is highly polar- 
ized in some cases (Zakamska et al.|2005l ). 15% ±5% qualify 
as radio loud (Zakamska et al. 12004 Lai fc Ho l2010l ). 

Because of the recent discovery, little is known about 
the possible exist ence of outflows in type 2 quasars. 
I Greene et alj l|201ll) analysed the kinematics of the extended 
ionized gas in a sample of 15 luminous SDSS type 2 quasars 
at z <0.5 in spatial scales of <10 kpc. They found kinemat- 
ically disturbed gas accross the host galaxies and proposed 
that AGN driven outflows are stirring up the gas on scales 
of kpcs. 

We investigate in this paper the presence of ionized gas 
outflows in a sample of 13 optically selected type 2 quasars at 
0.3< 2 <0.6. The observations and data analysis techniques 
are described in 52. The results are presented in 53 and 
discussed in 54. The results and conclusions are summarized 
in 55. 

We assume Q A =0.7, O M =0.3, H =71 km s" 1 Mpc" 1 . 



2 OBSERVATIONS AND ANALYSIS 

The data were obtained with the Focal Reducer and Low 
Dispersion Spectrograph (FORS2) for the Very Larg e Tele - 
scope (VLT) installed on UT1 (Appenzeller et al. 119981 ). 
The observations were performed in two different runs: 8 
and 9th Sept 2008 (4 objects) and 17 and 18 April 2009 (9 
objects). All spectra were obtained with the 600RI+19 grism 
and the GG435+81 order sorting filter. The useful spectral 
range was -5030-8250 A in the 2008 run and -5300-8600 
A in the 2009 run, so that in all cases at least the H/3 and 
[OIII]AA4959,5007 lines were within the observed spectral 
range. The pixel scales are 0.25" pixel -1 and 0.83 A pix -1 
in the spatial and spectral directions respectively. The spec- 
tral resolution, as measured from the sky emission lines, was 
FWHM=7.2±0.2 A and 5.4±0.2 A for the 2008 and 2009 
runs respectively. The slit width was 1.3" in 2008 and 1.0" 
in 2009. The observations, data reduction and object sam ple 
are described in detail in Villar-Martm et al. l|201lll2010l ). 

1-dim spectra were extracted from spatial apertures 
centered on the quasars spatial continuum centroid. The 
sizes were selected to optimize the S /N of the nuclear emis- 
sion lines, in particular for the detection of faint broad wings. 
In general, the aperture sizes are —1.5-2", corresponding to 
physical sizes —7-12 kpc depending on the object. Type 2 
quasars are often associated with extended emission line re- 
gions (EELR) which can exte nd sometimes for tens of kpc 
(e.g. Villar-Martm et al. 1201 ll ). well beyond the size of the 
NLR. The apertures we have defined contain emission from 
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both the NLR and, when existing, the EELR, although the 
contribution of the NLR is likely to dominate. 

The spectral profiles of [OIII] and H/3 were fitted with 
the STARLINK DIPSO package. One or more Gaussian 
components were used depending on the quality of the fits. 
This was evaluated based both on the errors of the fits and 
a visual inspection of the residuals. The [OIII] lines are 
in all cases the strongest and they were used to constrain 
the expected kinematic substructure of H/3. The ratio be- 
tween the [OIII]AA4959,5007 components was forced to be 
1:3, in accord with the theoretical value (Osterbrock 1989). 
Both components were forced to have the same FWHM and 
in general the separation in wavelength was also fixed. We 
found necessary to apply these constraints to all kinematic 
components in order to reject fits that were successful at re- 
producing the line shapes, but did not have physical mean- 
ing. 

To start, we attempted to fit the [OIII] profiles with 
a single Gaussian function, but this failed in all cases and 
very large residuals remained. The number of Gaussians was 
progressively increased until a good fit was achieved. Two 
or three Gaussians at most were required in all quasars. 
The H/3 line profile could be decomposed in different kine- 
matic components only in some objects where the line had 
enough signal and was not severely affected by atmospheric 
or galactic absorption. In those cases, it was assumed that 
the line consists of the same number of kinematic compo- 
nents as [OIII]. The FWHM and z of each component were 
constrained from those derived from [OIII] . Only those com- 
ponents with a detectable flux were considered in the final 
fit of H/3. Flux upper limits for the non detected components 
were also estimated. The FWHM (corrected for instrumen- 
tal broadening in quadrature), the velocity shift relative to 
the narrow [OIII] core (see (52.1), the flux, the flux relative 
to the total line flux and the [OIII] /H/3 ratio (when possible) 
were measured for all kinematics components. 

The kinematic measurements obtained for the first 5 
objects in Tables 1 and 2 are affected by slit effects (see 
Villar-Martm et al. 2011). The errors quoted take the re- 
sulting uncertainties into account. 

The results of the fits are shown in Table 2 and Figures 
1 to 13. 



2.1 The systemic velocity 

In order to characterize accurately the gas motions relative 
to their environment it is essential to determine the sys- 
temic redshift z sy3 of the galaxy. The safest way to do this 
is to measure the stellar redshift. This is often not possi- 
ble as for example in studies of type 1 AGNs, high redshift 
galaxies, etc. where the stellar continuum cannot be isolated 
from the bright AGN contribution and/or it is not detected 
with enough signal. In such cases, alternative methods have 
been applied. Under the assumption that the BLR gas fol- 
lows gravitational motions, the broad H/3 line centroid has 
been used to determine z sys . For type 2 objects, the [OIII] 
line fitted with a single Gaussian has often been used in- 
stead, both to measure z sys and the stellar velocity disper- 
sion. However, as different wo rks ha ve discuss ed (e. g. Greene 
et al . l201ll . I2009I . Husemann l201ll Boroson 12001 Boroson 
2005), this is risky, since the NLR kinematics often reveal 
non-gravitational motions (see SI). Indeed, the very asym- 
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Figure 1. Fit to H/3 and [OIII] AA4959.5007 for SDSS J2358-00. 
In each figure from 1 to 13 the data and the fit are shown in 
black and orange respectively. The same colour and line style 
are used for each kinematic component in all emission lines 
([OIII]AA4959,5007 and, when possible, also H/3). Red and blue 
are used for the most redshifted and blueshifted components re- 
spectively. Green is used for an intermediate redshift. When two 
components are identified with equal z within the errors (as the 
example in this figure), the same colour but different line style are 
used for each one. The vertical black line marks the wavelength 
of the very narrow [OIII] peak, which we have used to caculate 
the velocity shifts. 
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Figure 2. SDSS J0025-10. Line styles and colors as in Fig. 1. 



metric and sometimes broad [OIII] spectral profiles often 
make the determination of it s own z uncertain by up to 200 
km s _1 (e.g. Heckman et al. 119811 b 

There are three objects in our sample with enough sig- 
nal in the continuum and detected stellar features to try to 
estimate z sya : SDSS J0217-00, SDSS J1337-01 (both with 
strong stellar features) and SDSS J0123+00 (faint features). 
For this, we use the full spectral fit ting m ethod with the 
Starlight code (Cid Fernandes et al. l2005t ). several s ets of 
evolutionary stellar models (Gonzalez Delgado et al. 120051 ; 
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Figure 3. SDSS J0123+00. Line styles and colors as in Fig. 1. 
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Figure 5. SDSS J0234-07. Line styles and colors as in Fig. 1. 
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Figure 4. SDSS J0217-00. Line styles and colors as in Fig. 1. 



Chariot & Bruzual (in prep.)) and two of the most up to 
date high / inte rmedi ate spectral resolution stellar libraries 
(Martins et al. |2005l ; Sanchez-Blaquez et at. l2006t ). 

We find that in all three cases there is a variety of mod- 
els that produce acceptable fits (see one example in Fig. 14). 
These models allow a range of stellar redshifts which de- 
pend mainly on the single stellar population models used as 
templates for the fits and the weight given to the G band 
and Call K stellar absorption features. The allowed redshifts 
imply a velocity range of ~100 km s~\ which is the main 
uncertainty affecting the method. 

Although the [OIII] line in the nuclear spectra of our 
sample is in general quite broad and shows complex kine- 
matic substructure (see 83), it shows a very narrow core 
whose redshift can be estimated with an accuracy of several 
km s _1 . We will use this as a measurement of z aya . 

In the three objects mentioned above, it is found that 
this core is at <100 km s -1 of the stellar velocity, the un- 
certainty being due to the range of stellar redshifts allowed 



Figure 6. SDSS J0955+03. Line styles and colors as in Fig. 1. 
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Figure 7. SDSS J1153+03. Line styles and colors as in Fig. 1. 
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Figure 8. SDSS J1228+00. Line styles and colors as in Fig. 1. 
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Figure 11. SDSS J1407+02. Line styles and colors as in Fig. 1. 
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Figure 9. SDSS J1307-02. Line styles and colors as in Fig. 1. 
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Figure 12. SDSS J1413-01. Line styles and colors as in Fig. 1. 



SDSS1546-00 




6400 



6500 



6600 



6700 



680 6740 6800 6860 



6920 



Figure 10. SDSS J1337-01. Line styles and colors as in Fig. 1. 



Figure 13. SDSS J1546-00. Line styles and colors as in Fig. 1. 
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Figure 14. Fit to the continuum of SDSS J0217-00 using stellar 
population synthesis models. The spectrum has been translated 
to z=0. The fluxes have been normalized to the flux at 4020A. 
Different templates produce a good fit, so that a range of stel- 
lar (and therefore systemic) redshifts is allowed. These imply a 
velocity range of ~100 km s — 1 . 



stars in most quasar spectra studied here, we will use the 
narrow [OIII] core in all objects instead. Thi s method is 
supported by the results of I Greene fc Hd (|2005l ). They com- 
pared the gaseous and stellar kinematics in the nuclei of 
a large sample of active galaxies using spectra of some- 
what better resolution (R= A/AA ~1800 vs. R-900-1400 
at the [OIII] wavelength in our spectra, depending on the 
object). They found that the FWHM of the [OIII] core, i.e. 
after the asymmetric wings were removed, and FWHM„( ars 
are correlated (although with a large scatter) according to 
FWHM core =(1.24±0.76) x FWHM siars . 

Based on this result, they propose that FWHM core is a 
reasonable tracer of the stellar velocity dispersion. It seems 
reasonable to expect that it also provides a good estimation 
of t he stellar redshift an d therefore z sya . 

iGreene et al.l (|2009h questioned later the validity of this 
result for luminous obscured AGNs, based on the lack of cor- 
relation between the [OIII] FWHM and the stellar velocity 
dispersion. However, the analysis is based on single Gaussian 
fits to the emission lines that ignore their complex kinematic 
substructure. Given the high complexity of the line profiles 
(see 53.2), it is not surprising that no correlation is found. 
As Table 1 shows, the FWHM of the core is often noticeably 
smaller than that resulting from a single Gaussian fit to the 
line profile. 



2.2 The identification of kinematically perturbed 
gas 

We are specially interested in those gas components whose 
kinematic properties suggest non-gravitational, perturbed 
motions. In order to identify them, we need to compare 
the FWHM of each kinematic component with the stellar 
velocity dispersion FWHM s t ars . Similar values could be ex- 



plained if the gas is supported by random motion^, while 
rotation would produce narrower FWHM for the gas than for 
the stars. Broader gas components would imply kinematic 
perturbation. 

In general, it is not possible to measure FWHM„( ars 
directly from our spectra. Instead, following Greene & 
Ho (2005) we will infer the value for each object using 
FWHM core (52.1). Although the scatter of the correlation 
between these two parameters produces large uncertainties, 
the FWHM core values in Table 1 suggest FWHM sl „„ in the 
range ~250-900 km s _1 (Table 2), consistent with values 
m easure d in other type 2 quasars (Bian et al. |2007| . Liu et 
al. 120091 ') using stellar features. 

The velocity shifts Vs of the different [OIII] kinematic 
components relative to the narrow line core will provide ad- 
ditional information to characterize the kinematics. 



3 RESULTS 

3.1 Integrated line profiles 

We show in Table 1 the FWHM of [OIII] and H/3 lines. For an 
initial, general characterization of these quantities we have 
applied single Gaussian fits to both lines. The goal of this is 
to compare with similar studies of type 1 quasars. 

The [OIII] FWHM values in Heckman, Miley & Green 
( 1984) subsample of 35 type 1 quasars are in the range [260- 
850] km s _1 , with a median value of 460 km s _1 . If only the 
10 radio quiet quasars are considered the range is the same, 
but the median value increases to 570 km s -1 . The range in 
our sample is [410TT0 to 1860T20] km s _1 and the median 
is 510 km s" 1 . 

A visual inspection of the [OIII] spectral profiles re- 
veals t hat it is asymmetric i n the majority of objects. Fol- 
lowing [Heckman^et^al] l|l98ll ). we have quantified the degree 
of asymmetry at 10%, 20% and 50% of the line maximum 
intensity using the parameters AI10, AI20 and AI50. AI10 
is defined as ^g+^glg , where WL10 and WR10 are the 
half widths to the left and right of the center of the narrow 
[OIII] core at 10% intensity level. The widths were corrected 
for instrumental broadening. AI20 and AI50 are defined in 
the same way, but at 20% and 50% intensity levels respec- 
tively. AI10 is introduc ed co mpared with previous studies 
(see also Spoon & Holt [2010) because we find that in most 
objects the asymmetry is perceived at very low intensity lev- 
els, in the form of faint wings. The schematic representation 
of the various parameters involved in these calculations is 
shown in Fig. 15. The half widths WL10', WR10', etc have 
a " ' " because they correspond to values not corrected for 
instrumental broadening. 

AI10, AI20 and AI50 are shown in Table 1. Positive, 
negative and values indicate blue (highlighted in bold), 
red (italic fonts) and no asymmetries (normal fonts) respec- 
tively. It can be seen that the fraction of asymmetric profiles 
increases at lower intensity levels (7/13 at 50%; 9/13 at 20% 
and 10/13 at 10% levels). While at 50% level there is a vari- 
ety of asymmetries with no clear preference for a given sign, 

1 Type 2 quasars are most f requently associa ted w ith elliptical 
galaxies (Mainieri et al. l201ll Zakamska et al. |200(J) . so random 
motions are possible 
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430±10 


0.03±0.04 


0.00±0.03 


0.00±0.03 


SDSS 


J0025-10 


0.303 


435±55 


420±70 


320±80 


0.33±0.05 


0.36±0.08 


0.42±0.09 


SDSS 
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0.399 


445±45 


470±60 


430±50 


0.10±0.05 
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0.04±0.06 


SDSS 


J0217-00 


0.344 


910±30 


700±60 


700±60 


0.22±0.03 


0.06±0.04 
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SDSS 


J0234-07 


0.310 


510±60 


420±60 


450±50 


0.39±0.04 


0.52±0.05 


0.53±0.09 


SDSS 


J0955+03 


0.422 


1320±30 


N/A 


940±60 


0.13±0.03 


0.14±0.07 


-0.15±0.09 


SDSS 


J1153+03 


0.575 
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N/A 


440±20 
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-0.05±0.04 


-0.06±0.09 


SDSS 
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1860±20 


N/A 


HOOilOO 
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455±15 
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0.20±0.10 


SDSS 


J1413-01 


0.380 


550±10 


570±40 


500±10 


0.03±0.04 


0.03±0.03 


0.00±0.03 


SDSS 


J1546-01 


0.383 


500±40 
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430±20 
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0.43±0.08 


0.00±0.05 


Median 


0.383 
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0.13 


0.07 


0.00 
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440 


93 


230 


0.28 


0.24 


0.26 



Table 1. Type 2 quasar sample studied in this work with kinematic measurements based on single Gaussian fits to the [OIII] and H/3 
lines. The redshift is shown in column (2). The [OIII] and H/3 FWHM are are in columns (3) and (4). Column (5) gives the FWHM 
of the narrow [OIII] core that remains after the asymmetric line wings are removed. AI10, AI20 and AI50 in columns (6), (7) and (8) 
quantify the degree of asymmetry of the [OIII] line profile at 10%, 20% and 50% intensity levels respectively. Positive, negative and cero 
values (taking errors into account) correspond to blue, red and no asymmetries, which are highlighted in bold, italic and normal fonts 
respectively. 

cant fraction of the line emitting gas is moving radially in 
an outflow (see 81). 



3.2 Spectral decomposition of the H/3-[OIII] line 
profiles. 

The [OIII] spectral line profiles show complex kinematic sub- 
structure in all 13 objects: 3 kinematic components are re- 
quired to produce the best fit in 8/13 objects, while 2 are 
required for the remaining 5 (Fig. 1 to 13). 

All quasars in our sample have one or two kinematic 
components with [OIII] FWHM <FWHM st „ s (Table 2), 
i.e., they show no sign of kinematic disturbance. The pres- 
ence of two components which seem to follow gravitational 
motions is not surprising, since the emissions from the NLR 
(probably dominant) and possibly also the EELR are ex- 
pected within the spatial apertures used to extract the spec- 
tra. According to studies of narrow line radio galaxies, the 
EELR is expected to have different kinematic, physical and 
ionizatio n prop erties than the NLR (e .g. Tadhunter, Fosbury 
& Quinn ll989l . Robinson et al. 119871 ). 

All objects show, moreover, a very broad compo- 
nent with [OIII] F WHM>2 x FWHM stars . In general, this 
component is broader than the maximum possi ble value 
of FWHM star . s allowed by iGreene fc Hoi (|2005h correla- 
tion scatter (see 52.2.). Therefore, the uncertainty on 
FWHM S ( ars does not affect our conclusion. Thus, the NLR 
of all type 2 quasars in our sample contains highly perturbed 
gas, in addition to the more quiescent, ambient gas whose 
FWHM is consistent with gravitational motions. 

The FWHM of the perturbed gas components are in 
the range 780±30 to 2500±200 km s" 1 with a median value 
of 1280 km s" 1 . Their velocity shifts relative to the [OIII] 
narrow core Vs are in the range 10±40 and up to ~-620±140 
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Figure 15. Schematic representations of the various parameters 
used to quantify the degree of asymetry of the [OIII] line profile 
at differente intensity levels, following Heckman et al. (1981). 



at lower intensity blue asymmetries become predominant: at 
10% level 9/13 objects show a blue asymmetry, vs. 2 with 
red and 2 more which no asymmetry. 

Our results are consistent with 

iHeckman. Milev fc Greeiil j 19841) . who found a clear 
predominance of blue asymmetric [OIII] profiles in radio 
quiet and flat steep spectrum (compact) radio loud quasars 
at z ^0.5. The asymmetries disappear in steep spectrum 
(extended) radio loud objects. None of the objects in our 
sample are radio loud. 

An excess of blue asymmetries suggests that a signifi- 
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Figure 16. Velocity shift Vs in km s — 1 vs. R = F pyvHM~ )I11 ^ f° r 
all kinematic components isolated in the spectral fits. Coloured 
symbols correspond to the two peculiar objects suspected to have 
intermediate orientation between type 1 and type 2 (green for 
SDSS J1228+00; orange for SDSS J1153+03). Excluding these, 
there is a global trend for the perturbed gas components (R >2) 
to be blucsfhited, while more quiescent kinematic components 
R <1) show small or null velocity shifts. 



Figure 17. Velocity shift Vs (km s _1 ) vs. the FWHM (km s _1 ) 
in log of all kinematic components isolated in the the [OIII] line 
for all quasars in the sample, except the two peculiar objects 
SDSS J1228+00 and SDSS J1153+03. Equal symbols correspond 
to the same object. There is a general trend such that in each 
object the broadest component shows also the largest blucshift. 



km s , wit h a m edian value of ~-145 km s 1 (see also 
Greene et al. l201ll ). 

We have plotted R = F J^Jfl° t 111 ^ vs - Vs for all 
kinematic components (Fig. 16). The two objects high- 
lighted in colour are SDSS Jl 153+03 (oran ge) and SDSS 
J1228+00 (green 1 ). IvTUar-Martfn et al.l (|201ll ) proposed that 
these quasars might have an intermediate orientation be- 
tween type 1 and type 2, so that we have a direct view of an 
intermediate density region, more interior and closer to the 
black hole than the NLR we see in the rest of the sample. 
Fig. 16 shows that they are also peculiar from the kinematic 
point of view, since they are the only ones showing kinematic 
components with large red velocity shifts. 

If these two peculiar objects are excluded, Fig. 16 shows 
a strong trend for the highly perturbed gas (i.e. kinematic 
components with R >2) to be blueshifted, while the more 
quiescent gas (components with R <1.0) shows in general 
no or small velocity shifts. 

Vs is plotted vs. the FWFHM in Fig. 17 for the individ- 
ual kinematic components isolated in each quasar (the two 
peculiar quasars are excluded). The same symbol is used for 
each object. Selecting a given symbol, it can be seen that 
the largest FWHM values (i.e. the most perturbed gas in 
a given object) are in general associated with the largest 
blueshifts (see also Table 2). This result is found in all ob- 
jects where the velocity uncertainties allowed a definitive 
conclusion (7/11). 

Thus, the perturbed gas shows both broad FWHM com- 
pared with the stellar velocity dispersion and a trend to show 
the largest blueshift. The perturbed gas is therefore respon- 
sible for the blue line profile asymmetries. All together sug- 
gest the presence of ionized gas outflows in the inner <few 
kpc, such that the receiding gas is obscured by dust. 

The relative contribution of the perturbed (i.e. outfiow- 
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Figure 18. Relative contribution R2 to the total [OIII] line 
flux of the perturbed component (column 5 in Table 2) vs. 
R =FWHM/FWHM s ( ars for the perturbed kinematics compo- 
nents (i.e. the outflowing gas). An anticorrclation is observed. 
Kinematic components with more extreme kinematics (larger R) 
have a smaller relative contribution to the total line flux. 



ing) gas to the total [OIII] flux R2 — ^jjjpsr (column 
5 in Table 2) varies strongly from object to object. It can 
account for ~10% and up to ~70% of the total line flux. R2 
is plotted vs. R in Fig. 18. An anticorrelation is found such 
that, the more perturbed the outflowing gas is (larger R), 
the lower its relative contribution to the total [OIII] flux. 
A Spearman's rank correlation coefficient p—-0.79 is found 
and a statistical significance at the 99.95% level. 

The perturbed gas is also detected in H/3 in 7 objects. 
The [OIII]/H/3 ratio is in general very high (J?10, in most 
cases where it could be measured), which is similar or even 
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larger than that of the narrower kinematic components. This 
suggests very high excitation level of the outflowing gas. 

Given the size of the apertures used to extract the spec- 
tra, the effects of the outflows we see are constrained within 
radial spatial scales of r <few kpc from the central engine. 
Larger spatial scales of the tot al out flowing region are not 
discarded (e.g. Humphrey et al. I2010T ). 



4 DISCUSSION 

4.1 Outflows in the NLR of type 2 quasars. 

We have found evidence for outflows of ionized gas in the 
central regions (r <few kpc) of most type 2 quasars in 
our sample (11/13). Six objects in the sample were selected 
on the basis (among other things) of prior evidence of nu- 
clear kinematic perturbation (Villar-Martm et al. 2011) such 
as broad and/or highly asymmetric [OIII] lines. However, 
not obvious bias was applied to give preference to blue 
vs. red asymmetries. Moreover, six out of the seven ob- 
jects for which kinematic criteria were not applied also con- 
tain outflows. So, in spite of the apparent selection bias, 
the results suggest that the outflow phenomenon is com- 
mon in the NLR of type 2 quasa rs. This is consistent with 
iHeckman. Milev fc Greer] (jl984T ) . who found similar evi- 
dence in the NLR of the majority of r adio quiet type 1 
quasars (see also Leipski & Bennert I2006T ). 

The outflowing gas shows high ionization level with 
[OHl]/H/3 sometimes > 10, and often similar or larger than 
the more quiescent gas. A similar result has been found in 
some radio quiet type 1 quasars (e. g. Le ipski & Bennert 
l2006h and Seyfert galaxies (Veilleux Il99ll ). fndependently 
of the outflow origin, such high excitation suggests that 
stars are not responsible for ionizing the outflowing gas, but 
rather AGN related phenomena. This supports that the out- 
flows are within the quasar ionization cones, as commented 
above. Although contribution from the EELR is expected 
in the spectra of at least some objects, the line emission is 
likely to be dominated by the NLR so that the outflows are 
expected to be located in this region. 

The outflowing gas emits very broad lines 
(FWHM>2 x FWHM SIors , with FWHM in the range 
780±30 to 2500±200 km s" 1 ) and projected shifts relative 
to the systemic velocity in the range Vs — 10±40 and up to 
— -620±140 km s" 1 relative to the systemic velocity. The Vs 
are in general much less extreme than the motions implied 
by the large FWHM. Very broad lines can be produced by 
acceleration of clouds by or behind the bow s hock generated 
by the outflow (e.g. Villar-Martm et al. 1 19991). Also, it must 
be taken into account that the deprojected velocity shifts 
are likely to be larger, since these are type 2 objects and 
the cones axis is close to the plane of the sky. 

We have found a anticorrelation such that the strongest 
the kinematic perturbation (larger 7?), the lowest the rela- 
tive contribution of the outfowing gas to the total [OIII] line 
(smaller 722) . Since there is no obvious change in the ioniza- 
tion level of the outfowing gas relative to the ambient, non 
perturbed gas, this anticorrelation possibly reflects that, the 
more perturbed the gas is, the smaller the gas mass fraction 
the outflow involves. 

Based on the spatially resolved behaviour of the FWHM 



and Vs of [OI II] line in a s a mple of 15 type 2 quasars at 
0.11^ z $0.43. [Greene et al.l (|201 lh also proposed that out- 
flows of ionized gas are present in their sample. 



4.2 Outflow mass and mass outflow rate 

The mass of the ionized outflowing gas M can be estimated 
with: 

M _ 1A m p L (Hp) 
n e h vnp a e Jp 

where L(H/3) is the luminosity of the broad H/3 compo- 
nent, n e the density of the outflowing gas, m p the proton 
mass, a e Jp the effective H recombination coefficient for H/3 
(Osterbrock 1989) and hvufi is the energy of an H/3 photon. 

This refers to the outflowing gas detected within the 
slit. If the outflows are constrained to a region of size r <few 
kpc, we expect M to provide an adequate measurement of 
the t otal outflowing gas, otherwise (e.g. Humphrey et al. 
l20ld) the inferred masses are lower limits. 

Large uncertainties affect this calculation since the 
unreddened H/3 luminosity is unknown and we have no den- 
sity measurements. The spectral fits imply observed L(H/3) 
in the range (0.2-2.1) x 10 41 erg s" 1 . 

Using the [SII] double t from the SDSS spectra of type 
2 quasars, Greene et al. ()201 lh find a range of densities 
-250-500 cm" 3 , with a mean value of 335 cm" 3 . [OII]A3727 
is within the spectral range for SDSS J0217-00. The out- 
flow is also evident in this line, which presents a clear blue 
asymmetry. Given the clear detection of this wing, it is 
reasonable to think that n is below the critical density of 
this line n cr it=3xl0 3 cm" 3 . The gas densities estimated in 
other objects (radio galaxies, Seyferts, far infrared galax- 
ies with galactic superwinds) for the outflowing (nuclear or 
extended) gas cover a broad range from —several hundred 
cm" 3 (e.g Villar-Martm et al. 1 999, A rribas & Mediavilla 
1993, Heckman, Miley & Armus ll99oT ) to n e >10 4 5 (Axon 
et al. 1998). n =250 cm" 3 seems therefore a reasonable lower 
limit. Thus, we assume a range of densities n —250 - 3000 
cm" 3 for the outflowing gas. The implied values for M us- 
ing the uncorrected L(H/3) are in the range —(0.3-1.5) x 10 7 
Mq for n —250 cm" 3 and depending on th e object. The 
values are 120 times lower for n —3000 cm" 3 . iGreene et al.l 
l|201lf ) propose — 10 7 Mq as a strict lower limit assuming a 
gas density — fewxlOO cm" 3 . It must be taken into account 
that these authors argue that all the H/3 emitting gas is in- 
volved in the outflow (we find this is not the case), and that 
much lower densties are possible. 

For an illustrative comparison, the masses involved in 
neutral and molecular outflows in Seyfert and ultraluminous 
infrared galaxies (ULIRGs), which are constrained to spatial 
scales of <<few kpc, are > 10 8 M fll and often >10 9 M fl (e.g. 
Rupke, Veilleux & Sanders I2005I . Sturm et al. l201ll ). 

It is obvious that better constraints on the intrinsic 
L(H/3) and density of the outflowing gas are essential to 
obtain more tightly constrained masses. It is also essential 
to investigate the existence of outflows in other gas phases. 
Indeed, outflow studies in different types of galaxies suggest 
that the ionized phase reveals a very small fraction of the 
total outflowing mass, which is likely to be dominated by 
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Object 
(1) 


Asymm 
(2) 


Flux[OIII] 

X10- 16 
or g3 -l cm-2 
(3) 


[OIIl]i 


Hft t 


[01XI]i 


FWHM[OIII] 

km a" 1 
(7) 




[OIII]tot 


L[OIII] tot 
X 10 42 
erg s — 1 
(10) 


FWHM stars 

km S -l 
(11) 


\OIII\tot 

(4) 


UHtat 
(5) 


mh 

(6) 


V S 

km s- 1 
(8) 


HPtot 

(9) 


2358-00 


No 


75±6 
50±2 


0.60±0.06 
0.40±0.03 


0.62±0.0S 
0.38±0.05 


10±1.5 
11±1 


310±40 
1020±30 


10±40 
10±40 


10.9±0.4 


8.06 


350±40 



0025-10 Blue S.3±0.4 0.16±0.01 ^0.08 >9 <130 20±40 4.6±0.2 2.07 260±60 

35±1 0.67±0.03 ~1 3.4±0.1 500±40 -60±40 

8.8±0.9 0.17±0.02 <0.17 >4.3 1330±100 -80±45 



0123 + 00 


Blue 


5.9 + 0.9 
61.7+0.6 
8.2 + 0.7 


0.08 + 0.01 
0.81 + 0.02 
0.11 + 0.01 


<0.08 
~1 
^0.13 


>8.3 
8.8 + 0.2 
>7 


^200 
470 + 40 
1500+130 


20 + 50 
-15 + 50 
3 + 50 


10.3 + 0.3 


5.21 


350 + 40 


0217-00 


Blue 


30±2 
36 + 2 


0.46 + 0.03 
0.54 + 0.03 


0.63 + 0.09 
0.37+01.5 


10 + 2 
21 + 9 


470 + 50 
1530 + 40 


80 + 45 
-120 + 50 


14.6 + 2.3 


2.17 


570 + 50 


0234-07 


Blue 


32.5 + 0.9 
32±1 


0.48 + 0.02 
0.48 + 0.06 


0.50 + 0.04 
0.40 + 0.06 


11.1 + 0.9 

14 + 2 


310 + 70 
950 + 30 


+ 50 
-270 + 50 


12.6 + 0.6 


2.77 


360 + 40 


0955+03 


Blue 


0.9 + 0.2 
3.1 + 0.4 
2.6 + 0.3 


0.13 + 0.04 
0.47+0.08 
0.40 + 0.06 


Unk. 
Unk. 
Unk. 


Unk. 
Unk. 
Unk. 


470 + 80 
1110 + 80 
2500 + 200 


+ 20 
190 + 50 
-620 + 140 


16+4 


1.54 


740 + 50 


1153 + 03 


No 


26 + 1 
58 + 1 
9 + 2 


0.28 + 0.01 
0.62 + 0.02 
0.10 + 0.01 


Unk. 
Unk. 
Unk. 


Unk. 
14 + 4 
Unk. 


270 + 10 
730 + 40 
1450 + 80 


5 + 10 
-5 + 10 
330 + 70 


14.1 + 0.1 


15.7 


360 + 20 


1228+00 


Red 


6.4+0.3 
4.6 + 0.5 
22.2 + 0.4 


0.19 + 0.01 
0.14+0.2 
0.67+0.02 


Unk. 
Unk. 
Unk. 


Unk. 
Unk. 
Unk. 


550 + 20 
780+60 
2300 + 70 


-60+10 
640 + 30 
770 + 20 


12 + 1 


7.36 


890 + 80 


1307-02 


Blue 


17 + 2 
19 + 2 
13.6 + 0.7 


0.34+0.04 
0.40 + 0.04 
0.26 + 0.02 


0.23 + 0.03 
0.61 + 0.05 
0.16 + 0.04 


16 + 2 
7.0 + 0.8 
18 + 4 


250 + 10 
480 + 20 
1000 + 30 


-13 + 5 
56+9 
-190 + 20 


10.9 + 0.3 


3.21 


290 + 8 


1337-01 


Blue 


9.5 + 0.4 
22.3 + 0.5 


0.30 + 0.01 
0.70 + 0.02 


Unk. 
Unk. 


Unk. 
Unk. 


350 + 20 
1280 + 60 


30 + 6 
-290 + 10 


25 + 2* 


2.03 


440 + 40 


1407+02 


Blue 


55 + 1 
41 + 2 
15.0 + 0.6 


0.50 + 0.02 
0.37+0.02 
0.13 + 0.01 


0.57 + 0.03 
0.43 + 0.04 
<0.18 


9.4+0.5 
9.0 + 0.9 
>6 


290 + 10 
720 + 20 
1670+120 


10 + 3 
-100 + 10 
-70 + 30 


10.9 + 0.4 


2.03 


320 + 8 


1413-01 


No 


14.2 + 2 
45 + 2 
14+2 


0.20 + 0.02 
0.61 + 0.03 
0.19+0.02 


0.19 + 0.02 
0.56 + 0.04 
0.25 + 0.03 


9+1 
10.0 + 0.6 
7+1 


260 + 20 
600 + 20 
1190 + 50 


-20 + 10 
40 + 8 
-145 + 30 


9.6 + 0.4 


3.07 


400 + 8 


1546-01 


Blue 


7.6 + 0.1 
5.5 + 0.5 


0.58 + 0.02 
0.42 + 0.04 


0.71 + 0.08 
0.29 + 0.08 


6.5 + 0.7 
12 + 3.5 


340 + 20 
780 + 30 


00+20 
-270 + 35 


8.3 + 0.3 


5.21 


350 + 20 



Table 2. Results of the spectral decomposition analysis of [OIII] and H/3. (1): Object (shortened) name. (2): Sign of the [OIII] profile 
asymmetry implied by AI10 (Table 1). (3) Flux in units of xlO - 16 erg s~ 1 cm" 2 . (4) Fractional contribution of a given kinematic 
component to the total [OIII] flux. (5) Idem for H/3. (6) Ionization level as measured from the [OIIIJ/H/3 ratio. (7) FWHM. (8) Velocity 
shift Vs relative to the [OIII] narrow core. (9) [OIIIJ/H/3 calculated with the total line fluxes (rather than the individual kinematic 
components as in (6)). (10) Total [OIII] luminosity in units of 10 42 erg s _1 . (11) Stellar velocity dispersion calculated according to 
FWHM[OIII]core = (1.24±0.76)xFWHM at ars (Greene & Ho 2005). The errors do not account for the scatter of the correlation (see <53.2). 



the molecular and neutral phases. As an example, the mass 
of ionized gas in the (galactic) superwind region of the M82 
starburst galaxy is ~2xl0 Mq (Heckman, Armus & Miley 
1990), while the molecular mass involved in the outflow is 
>3xl0 8 Mq (Walter, Weiss & Scoville l2002h . 

The mas s outfl ow rate M can be calculated as (Dopita 
& Sutherland [2051 ): 

^ _ 3 v L(Hf3)m p _ 2.1 v M 
r n hv H pct e il r 

Assuming a spherical geometry, vo is the expanding ve- 
locity of the outflowing bubble and r its radius. In addition 
to those affecting M, additional uncertainties are involved 
for Vo, which is likely to be >V S (Table 2) and r. Although 
different works suggest that the radius o f the N LR in quasars 
is_~l-several kpc (e.g. Bennert et al. [20021 . Greene et al. 
l201ll ). this is only an upper limit for r, since the size of the 
outflowing region could be smaller. 

In spite of the uncetainties, it is interesting to study 
which conditions are required to produce different M. For 
instance, mass flow rates of the NLR outflows discovered in 
Seyfert and nearby narrow line radio galaxi es are typicially 
<few Mq yr _1 (Rupke, Veilleux & Sanders 2005, Heckman 



et al. 1981; see also Humphrey et al. 2010). Such values are 
certainly possible in our sample within a rather broad range 
of allowed parameter values. 

More restrictive conditions are required to produce 
M>50 Mq yr _1 as observed in some galactic scale outflows 
power ed by active star formation (e.g. Heckman, Ar mus fc 
Milev ll990l ; Lehnnert & Heckman [1999 ; Rupke et al. 120021 ). 
Let us consider a scenario which favours the efficiency of the 
outflow: low densities (n e =250 cm -3 ) and large line luminos- 
ity corrections (e.g. intrinsic L(H/^)^10xobserved L(H/3)). 
With these favourable prerequisites, the size of the outflow- 
ing region must be r <<lkpc. Studies of ionized gas outflows 
using UV and X ray absorption lines indicate that such small 
distances relative to he continuum source are p ossible. If the 
outflo w region has r ~several kpc as argued bv lGreene et al.l 
(2011), it is very unlikely that such high mass outflow rates 
are generated, since most of the outflowing gas should be 
collimated within a very narrow cone near the plane of the 
sky for all objects. 

It is clearly essential to determine the location of the 
outflows to constrain the mass outflow rates. 
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5 CONCLUSIONS 

We have analyzed the spatially integrated kinematic proper- 
ties of the NLR gas in a sample of 13 optically selected SDSS 
type 2 quasars at z ~0.3-0.6 using the [OIII]AA4959,5007 
lines. Optical long-slit spectroscopy at moderate spectral 
resolution R~900-1300 obtained with FORS2 on the VLT 
has been used for this purpose. 

The line profiles show a significant asymmetry in 11 of 
the 13 quasars. There is a clear preference for blue asymme- 
tries, which are found in 9/13 objects. An excess of blue line 
asymmetries have been observed in other active classes and 
non active galaxies. They are suggestive of outflows, such 
that the receding expanding gas is preferentially obscured. 
We propose the same interpretation. The outflows we see 
are constrained within spatial scales of r <few kpc from the 
central engine. 

This scenario is favoured by other results we find: in 
addition to more quiescent gas, whose kinematic properties 
are consistent with gravitational motions, we have discov- 
ered highly perturbed gas in all objects. This gas emits very 
broad lines (R — F ^ I I ^}° 111 ^ <^2). While the quiescent gas 
shows small or null velocity shifts relative to the systemic ve- 
locity, the highly perturbed gas trends to show larger shifts. 
These, moreover, are blueshifts in general. In each object 
the most perturbed gas trends to have the largest blueshift 
as well. All together support that the perturbed gas, which 
is responsible for the blue asymmetries on the line profiles, 
is outflowing. 

The projected velocity shifts relative to the narrow core 
of the [OIII] line, which we have considered a reasonable 
tracer of the systemic velocity, are in the range ~10±40 and 
up to ~-620±140 km s~\ with a median value of ~-145 km 
s _1 . The intrinsic expanding velocities are likely to be larger, 
given the orientation of the objects. We FWHM for the out- 
flowing gas in the range 780±30 to 2500±200 km s _1 with 
a median value of 1280 km s^ 1 . The relative contribution of 
the outflowing gas to the total [OIII] line flux varies from 
object to object in the range ~10-70%. An anticorrelation 
is found such that, the more perturbed the outflowing gas 
is (larger R), the lower its relative contribution is to the to- 
tal [OIII] flux . Our results suggest that outflows with more 
perturbed kinematics involve a smaller fraction of the total 
mass of ionized gas. 

The outflowing gas is in general very highly ionized, as 
suggested by the large [OIII]/H/3 values (often >10). Inde- 
pendently of the outflows origin, this suggests that the gas 
is ionized by the quasar, rather than stars, an therefore it is 
located within the quasar ionization cones. 

We argue that our results suggest that nuclear dusty 
outflows are a common phenomenon in optically selected 
type 2 quasars at 0.3< z <0.6. 
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